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The ability to accurately model details of inlet back flow for inducers operating at low- 
flow, off-design conditions is evaluated. A sub-scale version of a three-bladed liquid 
hydrogen inducer tested in water with detailed velocity and pressure measurements is used 
as a numerical test bed. Under low-flow, off-design conditions the length of the separation 
zone as well as the swirl velocity magnitude was under predicted with a standard k-e model. 
When the turbulent viscosity coefficient was reduced good comparison was obtained at all 
the flow conditions examined with both the magnitude and shape of the profile matching 
well with the experimental data taken half a diameter upstream of the leading edge. The 
velocity profiles and incidence angles at the leading edge itself were less sensitive to the back 
flow length predictions indicating that single-phase performance predictions may be well 
predicted even if the details of flow separation modeled are incorrect. However, for 
cavitating flow situations the prediction of the correct swirl in the back flow and the 
pressure depression in the core becomes critical since it leads to vapor formation. The 
simulations have been performed using the CRUNCH CFD® code that has a generalized 
multi-element unstructured framework and an advanced multi-phase formulation for 
cryogenic fluids. The framework has been validated rigorously for predictions of 
temperature and pressure depression in cryogenic fluid cavities and has also been shown to 
predict the cavitation breakdown point for inducers at design conditions. 

Nomenclature 

C m , c g , C L Speed of sound in mixture, gas, liquid 

v Viscous flux vector 

F,G Flux vectors 

„ Mixture enthalpy 

B Turbulent kinetic energy, turbulent dissipation rate 

y ,K ^ Rate constants for vapor and liquid formation 

\ Rate of vapor mass formation 

Pressure 

m'PgiPL Density of mixture, gas, liquid 

, E, S Vector of conservative variables 
v Vector of primitive variables 

Source terms 

S g Cavitation source term 

T sat Temperature saturation 

u,v,w x,y,z - components of velocity 

D H v Latent heat formation of vapor 

<p g (/) L Void fraction of gas, liquid 

&& Ci and C 2 Modeling constants 
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I. Introduction 

O UR interest is in the development of a computational framework to simulate cavitating liquid rocket 
turbomachinery that employ cryogenic working fluids. Liquid rocket systems are a subset of a broader class of 
pumps (e.g. refrigerant systems, boiler feed pumps, etc) where the operating temperature is elevated relative to the 
critical temperature of the fluid and thermodynamic effects of cavitation play an important role. At these operating 
temperatures, the ratio of liquid to vapor density is lower and consequently more liquid mass has to vaporize to 
sustain a cavity. Therefore evaporative cooling effects are more pronounced and result in the lowering of the mean 
fluid temperature in the cavitating region. Since the fluid thermodynamics properties (i.e. vapor pressure, density) 
are a strong function of temperature at these conditions, thermal effects suppress cavitation and lower the cavity 
pressure in a mean sense. Typically this results in improved mean performance of cryogenic pumps; liquid 
hydrogen systems being an extreme example where the pump may continue to generate head even when the fluid is 
boiling at the inlet. 

The thermal effects of cavitation were studied extensively by numerous researchers through the 1970’s 
including: Hord \ Stahl and Stepanoff 2 , Ruggeri and Moore 3 , Holl 4 , and Brennen 5 among others. Stahl and 
Stepanoff 2 were the first to estimate head depression (A// v ) values due to thermodynamic effects using the so- 
called ‘B-factor’ method based on a quasi-static theory where the temperature depression was estimated in terms of 
the ratio of the vapor volume to liquid volume. They provided a graph to evaluate NPSH corrections for 
hydrocarbons based on this methodology. More elaborate correlations, which included dynamic effects were given 
later by Hord \ Ruggeri and Moore 3 , and Holl 4 . They collected extensive experimental data of cavity pressure and 
temperature depressions for a variety of model shapes and fluids and correlated the results using variants of the B- 
factor theory. The semi-empirical procedures outlined by Ruggeri and Moore 3 continue to be used as an 
engineering tool for predicting the thermodynamic depression in pumps. A more rigorous numerical procedure was 
developed by Cooper 6 where a baratropic equation of state was used to define the two-phase mixture and thermal 
effects were evaluated with a resulting non-dimensional vaporization parameter. Most of these techniques, however, 
require some degree of empiricism. Therefore, from a more fundamental modeling perspective, this discussion 
highlights the need for a generalized formulation that takes the energy balance into account when simulating 
cavitation for cryogenic flows. 

For turbomachinery simulations, the majority of cavitating pump simulations presented in the literature are 
limited to simulating idealized liquids (no thermodynamic effects) at design conditions where large-scale inlet back- 
flow is not present. Typical simulations show comparisons with data at design conditions for the head coefficient 
and the critical N ss number at which performance breakdown occurs (Hosangadi et al. 7 , Athavale and Singhal 8 , 
Dupont and Okamura 9 , Medvitz 10 ). However flows at off-design conditions, where large scale unsteadiness and 
high dynamic pressure loads are observed, cannot at this point be reliably predicted Simulation of cavitation 
instabilities and rotational cavitation modes in pumps have not been simulated by any group to the best of our 
knowledge. 

The development of a framework for cavitating, cryogenic fluids with real fluid property variations was 
addressed in our earlier work (Hosangadi and Ahuja n ) and is an extension of our model for idealized fluids (Ahuja 
et al. 12 ). To simulate this class of flows, a generalized multi-phase formulation has been developed that rigorously 
models thermal effects of phase change and the accompanying property variations. Thermal equilibrium is assumed 
and fluid thermodynamic properties are specified along the saturation line using the NIST-12 databank. The multi- 
phase model has been implemented in the CRUNCH CFD® code, which is a multi-element based unstructured code 
’ . The underlying philosophy in the CRUNCH CFD® code is to tailor the grid topology to resolving the dominant 
flow phenomena and the structural complexity of the problem. This is achieved by utilizing a combination of 
hexahedral, tetrahedral, prismatic and pyramidal elements in mesh construction. Such a framework is particularly 
attractive for complex turbomachine configurations, since high quality grids can be generated very efficiently with 
minimum skewness. 

The multi-phase, real-fluid formulation was rigorously validated for cavitation in liquid nitrogen and hydrogen 
(Hosangadi and Ahuja n ). Predictions of temperature and pressure depression for flow over hydrofoils were 
compared with detailed experimental data obtained by Hord 1. The framework was subsequently applied to the 
study of the full-scale LOX inducer (Hosangadi et al. 15 ) whose cavitating performance at design conditions was 
compared to that of sub-scale water test. The increase in suction performance of the LOX pump resulting from 
thermal effects was identified at design flow conditions where large scale flow separation is not present and the flow 
is relatively steady. 

In the present effort, our interest is in the simulation of flowfields at off-design, low-flow conditions where 
large scale separation is present and a swirling back flow extends upstream of the leading edge. The backflow 
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results in the acceleration of the meridional velocity in the core and the generation of a radial pressure gradient with 
a low pressure core to sustain the swirl. Under cavitating conditions the separated flow can interact with the phase 
change mechanism to generate cavitation surge instability wherein large-amplitude axial fluctuations of pressure and 
mass flow thru the pump may be observed. While the simulation of cavitating flow at these off-design conditions is 
our ultimate goal, the first step is to validate the ability to predict the single-phase back flow. 

Here we focus on evaluating the ability of a Navier-Stokes methodology with a two-equation k-e model to 
accurately predict the details of the single-phase flowfield upstream of the leading edge. The inducer chosen for our 
study was a three-bladed liquid hydrogen inducer designed at NASA Marshall (Thornton 16 ) using the Agile 
Engineering Design System tool of ConceptsNREC, and subsequently tested at ConceptsNREC in water for a 
subscale configuration (Japikse and Baun, 17 ). Detailed measurements of velocity and pressure are available 
upstream of the leading edge while the flow deviation angle downstream of the trailing edge was also measured. As 
we shall discuss later when backflow was present, the standard k-e turbulence model overpredicted the turbulence 
viscosity in the separated region and hence underpredicted the back flow region. Good comparison was obtained 
when the turbulence constants were altered to reduce the level of turbulent viscosity indicating the sensitivity of the 
separation zone to this quantity. We are currently involved in evaluating if other turbulence models such as the 
realizable k-e model perform better under these flow conditions. 

II. Multi-Phase Equation System 

The multiphase equation system is written in vector form as: 


dQ dE dF dG 
dt dx dy dz 


= S + D V 


( 1 ) 


The multiphase equation system is written in vector form as: 

Here Q is the vector of dependent variables, E, F and G are the flux vectors, S the source terms and D v 
represents the viscous fluxes. The viscous fluxes are given by the standard full compressible form of Navier Stokes 
equations. The vectors Q, E and S are given below with a detailed discussion on the details of the cavitation source 
terms to follow later: 
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Here, p m and h m are the mixture density and enthalpy respectively, and (|> g is the volume fraction or porosity of the 

vapor phase. The mixture energy equation has been formulated with the assumption that the contribution of the 
pressure work on the mixture energy is negligible which is a reasonable assumption for this flow regime. The 
source term for the vapor phase arises from rate of vapor mass generation due to cavitation m t and the 
corresponding source term for the energy equation is given as m l • h fg where h fg is the change in enthalpy resulting 

from the phase change and is a function of the local fluid temperature. 

The mixture density and gas porosity are related by the following relations locally in a given cell volume: 


= r gf g + r if l 

(3) 

1 = <Pg + <I>L 

(4) 
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where p g , p L are the physical material densities of the gas and liquid phase respectively and in general are functions 
of both the local temperature and pressure. 

Thus far we have not made any statements defining the temperatures characterizing the liquid and vapor. In 
general, the liquid and vapor may not be in equilibrium locally and can have independent temperatures. Examination 
of temperature and pressure data for cavitation in Freon by Ruggeri 18 reveals that the saturation vapor pressure 
corresponding to the local fluid temperature in fact matches the local pressure measurement. This indicates local 
thermodynamic equilibrium that is exploited to significantly simplify Eqn. (1). The thermodynamic properties of 
the liquid and vapor in the cavity may now be defined by a single variable; the saturation temperature T sal . Hence, 
all thermodynamic properties (density, vapor pressure, viscosity, etc) of both the liquid and the vapor phase may be 
generated as a tabular function of the saturation temperature. In our study here, these properties were generated 
from the Standard thermodynamic database 12 available from NIST for pure fluids. 

The equation system as formulated in Eqn. (1) is very stiff since the variations in density are much smaller than 
the corresponding changes in pressure. Therefore to devise an efficient numerical procedure we wish to transform 
Eqn. (1) to a pressure based form where pressure rather than density is the variable solved for. An acoustically 
accurate two-phase form of Eqn. (1) is first derived, followed by a second step of time-scaling or preconditioning to 
obtain a well-conditioned system. We begin by defining the acoustic form of density differential for the individual 
gas and liquid phase as follows: 


dp=XrdP,dp^\dP (5) 

C g C L 

Here c„ is the isothermal speed of sound fj?Lj in the pure gas phase, and c L is the corresponding isothermal speed 

\dPgJj 

of sound in the liquid phase, which is a finite-value. We note that in Eqn. (5) the variation of the density with 
temperature has been neglected in the differential form. This assumption was motivated by the fact that the 
temperature changes are primarily due to the source term and not by the pressure work on the fluid i.e. the energy 
equation is a scalar equation. This simplifies the matrix algebra for the upwind flux formulation significantly, at the 
potential expense of numerical stability in a time-marching procedure. However, more importantly, there is no 
impact on the accuracy since the fluid properties themselves are taken directly from the thermodynamic data bank 
for each fluid. 


Following the discussion above, the differential form of the mixture density p m using Eqn. (5) is written as, 


dp m =(p g -p L ) d< t> g +\dP 

C 4> 

2 „ 2 r „ 2 


(6) 


Here, c 0 is a variable defined for convenience and is not the acoustic speed, c m , in the mixture, which will be defined 
later. Using Eqn. (6), Eqn. (1) may be rewritten as: 
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The numerical characteristics of the Eqn. (7) are studied by obtaining the eigenvalues of the matrix, 
The eigenvalues of the system are derived to be: 


dQ, 


A = (u + c m , u-c m , u,u,u,u,u) (9) 

where c m turns out to be the well-known, harmonic expression for the speed of sound in a two-phase mixture and is 
given as: 


1 

2 ~Pm 
c m 

The behavior of the two-phase speed of sound indicates that at either limit the pure single-phase acoustic speed 
is recovered. However, away from the single-phase limits, the acoustic speed rapidly drops below either limit value 
and remains at the low-level in most of the mixture regime. As a consequence, the local Mach number in the 
interface region can be large even in low speed flows. To obtain an efficient time-marching numerical scheme, 
preconditioning is now applied to the system in Eqn. (7), in order to rescale the eigenvalues of the system so that the 
acoustic speeds are of the same order of magnitude as the local convective velocities. 

The effects of turbulent mixing are accounted for by employing a two-equation k-e model. The turbulent 
viscosity is obtained by solving transport equations for the turbulent kinetic energy and its dissipation rate as 
follows, 
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where, G k , g& Cj and C 2 are modeling constants given as follows: 

c M =0.09, (T k =1.4, a e =1.4, 

Cj = 1.44, c 2 =1.92, c 3 = 2.9556 


( 12 ) 


For the inducer calculations presented here a wall-function procedure was used to specify the boundary 
conditions for the turbulence kinetic energy and the dissipation rate at the blade surfaces. 

A. Cavitation Source Terms 

In the present effort, the cavitation source term is defined via a simplified non-equilibrium, finite rate form as 
follows: 


m t = K /Pl0l + K bPA 


( 13 ) 


where the constant K f is the rate constant for vapor being generated from liquid in a region where the local pressure 
is less than the vapor pressure. Conversely, K b is the rate constant for reconversion of vapor back to liquid in 
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regions where the pressure exceeds the vapor pressure. Here, the rate constants are specified using the form given 
by Merkle 19. 
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We note that for steady attached cavitation this simplified form may be adequate since the cavitation time scales 
do not interact with the fluid time scales if the cavitation rate constants are fast enough. For unsteady cavitation 
modeling, however, it becomes essential to integrate bubble dynamics within a dense cloud framework wherein both 
the number density and mean local radius of the bubbles in an evolving cloud are tracked. The development of a 
more rigorous non-equilibrium source term model is a topic of ongoing research. 

III. Validation Study for Cavitation in Cryogenic Fluids 

The multi-phase formulation described above has been validated extensively by simulating experiments by 
Hord 1 on a cavitating hydrofoil. Hord performed sub-scale tests using both liquid nitrogen and hydrogen in a blow- 
down tunnel. The details of the tunnel and the hydrofoil geometry are given in Figure 1. The tunnel width is 1 inch 
while the hydrofoil width is 0.312 inches. 


HY W-CI; 

S *«**'<>« 



Quarter Caliber Hydrofoil with a Tapered Diffuser 



Tunnel Geometry 

Figure 1. Details of Experimental Set-up for Cavitating Hydrofoil (Taken from Hord, 1). 


Hence considerable blockage effects from the tunnel wall 
are present and this necessitated modeling the tunnel 
geometry in the simulations. To ensure that the tunnel 
blockage effects were being correctly modeled we performed 
single-phase, non-cavitating simulations and have compared 
it to Hord’s non-cavitating data in Figure 2. Excellent 
agreement is obtained giving confidence that the tunnel 
interaction is being captured. We note that the single-phase 
solution is insensitive to the Reynolds number since both 
liquid hydrogen and nitrogen show identical pressure 
profiles. 



Figure 2. Non-Cavitating Pressure 
Distribution on Hydrofoil with Tunnel 
Blockage Modeled (Hord 1 ). 
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Cavitating results in liquid nitrogen for one flow condition are discussed here. We refer the reader to Hosangadi 
and Ahuja 11 for additional validation studies in both liquid nitrogen and hydrogen. The operating range of liquid 
nitrogen varies from roughly 70K -100K. The variation of properties along the saturation line for liquid and vapor 
densities as well as the vapor pressure are shown in Figure 3. At 89 K, the slope of vapor pressure curve indicates a 
16 KPa increase for a 1 K change in temperature. 

The general characteristics of the cavitating flowfield are shown in Figure 4(a-d) for tunnel conditions in Run 
290C (Table I). Figure 4a shows the temperature profile in the cavity. The strong temperature depression at the 
leading edge of the cavity is evident with the gradual temperature recovery due to condensation in the rear of the 
cavity. The vapor volume fraction is qualitatively compared with a typical flow visualization of the flow (note that 
the flow conditions at which this visualization was done are not reported in the Hord’s report). The overall shape 
and features of the cavity appear to be similar. The computed pressure field indicates strong interaction between the 
cavity and the tunnel wall which is expected due to the relative scales of the geometry. 


Table I. Run Conditions For Liquid Nitrogen Cases 


Run 

Number 

Freestream 
Temp. (K) 

Freestream 
Vel. (M/S) 

Cavitation 

Number 

Cavity 

Length (CM) 

289C 

88.64 

23.5 

1.55 

2.29 

290C 

83.06 

23.9 

1.70 

1.9 

293A 

77.64 

24.0 

1.75 

1.52 

294F 

77.94 

9.8 

1.78 

1.52 




Figure 4. Cavitating Hydrofoil Flowfield For conditions of Run 290C in Liquid Nitrogen. 
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The quantitative comparisons of pressure and temperature depression in the cavity are compared with 
experimental data in Figure 5. Note that the pressure values plotted are (P - P v oo ). For a non-cryogenic case this 

value would be zero in the cavity, while values below zero in the cryogenic case indicate pressure depression due to 
thermal effects. In general excellent comparison is obtained for the leading edge temperature depression of 
approximately 2.5 K. The temperature recovery within the cavity compares well both in their slope as well as the 
length of the cavity. The temperature rise in the cavity closure region shows some differences; the computed 
solution recovers to the freestream value more quickly than does the data, which doesn’t quite fully recover to the 
freestream value. As per the discussion by Hord \ this probably was due to the unsteady effects in the cavity closure 
region whereby the thermocouples were not always enclosed in vapor giving erroneous readings. The 
instrumentation error given for the chromel-gold thermocouples is 0.20 K and the computed results are within the 
uncertainty of the experiments. 



Figure 5. Comparison of Cavitating Hydrofoil Simulations with Data for Liquid Nitrogen. 


The comparison of the pressure depression in Fig. 5a also indicates excellent overall comparison with data and 
within the instrumentation error bar of 0.69 N / cm 2 . The leading edge pressure depression is 23.5 percent relative 
to the freestream vapor pressure. It further illustrates why temperature effects have a substantial impact on the 
performance of cryogenic pumps. The experimental data plotted includes both the actual pressure measured 
(symbol: circle) as well as the saturation pressure values (symbol: square) corresponding to the temperature 
measurements. The close match between the actual pressure and saturation pressure values indicates that the 
thermodynamic equilibrium assumption is valid for liquid nitrogen flows. . 

IV. Inducer Simulations at Low-Flow Conditions 

The inducer chosen for our validation study is a three-bladed liquid hydrogen inducer designed at NASA 
Marshall (Thornton 16) using the Agile Engineering Design System tool developed by ConceptsNREC. The 
subscale configuration (0.37 scale) tested in water b ConceptsNREC had a diameter of 4.394 inches, a tip clearance 
of 0.014 inches, a shaft speed of 3840 rpm, and a flow rate of 195 gpm at the design point. Figure 6 illustrates the 
inducer blade angle, inducer wrap and blade thickness distribution as a function of meridional length. 

Figure 7 shows the three-dimensional geometry along with features of the multi-element unstructured grid that 
was used. In general, hexahedral cells are used around each blade while the region in between the hexahedral blocks 
is filled with tetrahedral cells. This strategy yields high-quality grids with minimal skewness even for extremely 
shallow blade angles that are typically associated with inducer blades. 

We present simulation results for three flow rates; 235 gpm (120% of design), 195 gpm (design), and 155 gpm 
(80% of design). Note that the inducer test showed large scale separation even at the design condition of 195 gpm 
while the 80% case obviously had a much larger separation zone. No backflow is observed at the 120% case of 235 
gpm. Detailed radial profiles of both the meridional velocity and swirl velocity as well as the pressure were 
measured 2.44 inches upstream of the leading edge tip location. The upstream experimental data was used to 
estimate the velocity profiles at the leading edge of the blades by correcting for the hub area ad by preserving 
conservation of mass and angular momentum along stream tubes (Thornton l6 ). Note that in the instances where we 
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compare simulations with the leading edge “data” the simulation profiles are the actual computed solution, while the 
“data” is in fact the estimate from the raw data half a diameter upstream. 
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Figure 6. Inducer Blade Design Parameter Variation with Meridional Length: a) Blade angle; b) Blade 
Wrap; and c) Blade Thickness (Taken From Ref. 16). 



Figure 7. Three-Dimensional Inducer Configuration and Multi-Element grid Topology. 

We begin by presenting results for the 235 gpm (120%) case for which no back flow is observed. The velocity 
and pressure profiles upstream in the inlet are uniform and hence we compare the profiles for the meridional 
velocity and incidence angles at the leading edge with the approximate 1-D procedure in Figure 8. Good 
comparison is obtained for both quantities and as we shall see from comparisons with the upstream data for other 


9 

American Institute of Aeronautics and Astronautics 



flow conditions, the minor differences observed are more likely due to the 1-D approximations used to process the 
experimental data. 



Meridional Velocity (ft/s) incidence (Degrees) 

(a) (b) 

Figure 8. Blade LeadingEdge Velocity Profiles for 120% of Design Flow Rate: 
a) Merdional Velocity; b) Incidence Angle. 


The results for the design condition of 195 gpm indicate that separation occurs at the leading edge with a 
substantial back flow region as is evident from the flow contours of the meridional velocity, pressure and the 
effective turbulent viscosity shown in Figure 9. A comparison of the radial flow profiles in Figure 10 with 
experimental data indicate that while the qualitative nature of the profiles are similar the extent of the back flow is 
underpredicted resulting in lower swirl velocities and higher values for pressure in the core relative to the shroud. 
The smaller extent of the back flow was found to result from an over prediction of the turbulent viscosity from the 
turbulence models that damped the reverse flow and not from numerical issues such as grid resolution and location 
of inflow boundary. 


Meridional Swirl 



Figure 9. Flow Contours for 100% of Design Flow Rate; a) Meridional Velocity, b) Swirl Velocity, c) 

Pressure, d) Turbulent Viscosity. 

A sensitivity study to the turbulent viscosity levels was undertaken by reducing the constant . As the constant 

is reduced from 0.09 to 0.06 and then to 0.03 we observe (from Figure 10) that the meridional velocity and swirl 
velocity profiles compare well with the experimental data while the pressure comparison is reasonable. The flow 
contours for c M of 0.03 indicate a separation zone almost twice the original separation length and turbulent 
viscosity levels that are roughly half the baseline case (Figure 1 1). 
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The meridional velocity and incidence angle profile at the leading edge are plotted in Figure 12 and they reveal 
an extremely interesting point: under non-cavitating conditions the velocity profiles and incidence angles at the 
blade leading edge are relatively insensitive to errors in the prediction of the separation zone. These results also 
indicate that the 1-D stream tube approximations to obtain the profile at the leading breakdown breakdown when 
separation occurs. The insensitivity of the leading edge profiles to the separation zone length would indicate that the 
single-phase performance predictions would not be affected significantly by errors in the estimation of the back flow 
but clearly cavitation performance would be quite sensitive to these errors since the extent of separation impacts the 
pressure depression in the core. Studies with other turbulence models that may perform better under swirling 
conditions (e.g. realizable k-e model) are currently underway. 

Simulations at 155 gpm (80% of design) were carried out with the c M value at 0.03 which yielded the best 

results at the design flow rate. Figure 13 shows the corresponding flow contours for meridional velocity, swirl 
velocity, and pressure. The back flow zone is substantially larger and extends nearly 3 diameters downstream. 



Figure 10. Radial Flow profiles compared with Data for 100 % of Design flow Rate Half a Diameter 
Upstream of Leading Edge: a) Meridional Velocity; b) Swirl Velocity; and c) Pressure. 


Meridional 
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(a). Cp=0.09 (b). Cp=0.03 

Figure 11. Sensitivity of Flow Contours to Turbulent Viscosity a): Cp=0.09; b) Cp=0.03. 
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Figure 12. Blade Leading Edge profiles for 100 % of Design Flow Rate: a) Meridional Velocity; b) Incidence 

Angle. 

The radial profiles for the velocity and pressure are compared with experimental data in Figure 14. Excellent 
comparison is obtained for the pressure and swirl components in particular. The back flow has high swirl near the 
shroud of approximately 18 ft/s which in turn generates a pressure gradient of 2.5 psi between the shroud and the 
core. This low pressure core will clearly result in a cavitating core as the inlet pressure is dropped and possibly lead 
to cavitation induced axial instabilities. Finally we present the head coefficient computed for the three flow rates 
modeled with the experimental data in Figure 15. Good comparison is obtained for all three cases indicating that the 
inducer performance is being simulated well. 



Figure 13. Flow Contours for 80% of Design Flow Rate: a) Meridional Velocity; b) Swirl Velocity; 

c) Pressure; d) Turbulent Viscosity. 



Figure 14. Radial Flow profiles compared with Data for 80% of Design flow Rate Half a Diameter Upstream 
of Leading Edge: a) Meridional Velocity; b) Swirl Velocity; and c) Pressure. 
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Figure 15. Computed Head Coefficient Compared with Experimental Data for Various Flow Coefficients. 

V. Conclusion 

A generalized multi-phase formulation has been developed for cavitation in “real” fluid flows (e.g., cryogenic 
fluids) that operate at temperatures close to their critical temperatures. It has also been previously shown to predict 
the cavitation breakdown point for inducers at design conditions where substantial back flow is not present. 

The focus of the current effort is in continued validation of this tool for predicting off-design performance of 
inducers. A sub-scale version of a three-bladed liquid hydrogen inducer tested in water with detailed velocity and 
pressure measurements was used as a numbered test bed. Under low-flow, off-design conditions it was found that 
the length of the separation zone as well as the swirl velocity magnitude was under predicted with a standard k-e 
model and the turbulent viscosity constant had to be reduced. With the modified constant good comparison was 
obtained at all the flow conditions examined with both the magnitude and shape of the profile matching well with 
the experimental data taken half a diameter upstream of the leading edge. The velocity profiles and incidence 
angles at the leading edge itself were less sensitive to the back flow length predictions indicating that single-phase 
performance predictions may be well predicted even if the details of flow separation modeled are incorrect. 
However, for cavitating flow situations the prediction of the correct swirl in the back flow and the pressure 
depression in the core becomes critical since it leads to vapor formation in the core that can couple with the 
hydrodynamics and generate cavitation surge. 
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